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ABSTRACT 

We explore the role of anisotropic thermal conduction on the evolution of supernova 
remnants through interstellar media with a range of densities via numerical simula- 
tions. We find that a remnant expanding in a dense environment can produce centre- 
bright hard x-ray emission within 20 kyr, and centre-bright soft x-ray emission within 
60 kyr of the supernova event. In a more tenuous environment, the appearance of a 
centre-bright structure in hard x-rays is delayed until about 60 kyr. The soft x-ray 
emission from such a remnant may not become centre bright during its observable 
lifetime. This can explain the observations that show that mixed- morphology super- 
nova remnants preferentially occur close to denser, molecular environments. Remnants 
expanding into denser environments tend to be smaller, making it easier for thermal 
conduction to make larger changes in the temperatures of their hot gas bubbles. We 
show that the lower temperatures make it very favorable to use high-stage ions as 
diagnostics of the hot gas bubbles in SNRs. In particular, the distribution of O Vlll 
transitions from shell-bright at early epochs to centre-bright at later epochs in the evo- 
lution of an SNR expanding in a dense ISM when the physics of thermal conduction 
is included. 
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1 INTRODUCTION 

Supernovae play an important role in driving the ther- 
mal evolution of the interstellar medium (ISM henceforth). 
The interiors of supernova remnants (SNR henceforth) are 
filled with extremely hot plasma, with temperatures of 
10^ — 10^ K. The classical picture of the expansion of a 
supernova remnant into the interstellar me dium identifie s 
four stages that are still germane today f see IWoltieJIlQ?^ ) . 
At very early times, the remnant is freely expanding into the 
undisturbed medium. After the shock sweeps up more mass 
than was ejected by the sup ernova, the r emnant expands 
adiabatically according to the lSedovl (|l959l) solution. When 
the cooling time of the matter behind the shock becomes 
comparable to the dynamical time of the expansion, a cold, 
dense shell is formed. Eventually, the expansion velocity of 
the shock becomes comparable to the turbulent velocity in 
the ISM, and the remnant disperses. 

SNRs can be classified into s everal morphologica l types; 
see Weiler & Sramek (lOSl) and Ijones et al.l (|l998l) for ex- 
cellent reviews. The classical or shell-type SNRs are the most 
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common and are shell-bright in both x-rays and radio emis- 
sion. The former arises from thermal bremsstrahlung and 
inner-shell metal line emission in the hot gas behind the 
shock. Synchrotron emission from cosmic rays and magnetic 
field swept up by the shock leads to the radio morphology. 
Plerionic remnants are powered by a central pulsar and have 
bright, nonthermal radio emission in their interior. Finally, 
there is a composite class, including radio and x-ray com- 
posite morphologies, which exhibit both a radio-bright shell 
and centre-bright x-ray emission. The x-ray composite rem- 
nants (also known as thermally composite or mixed mor- 
phology, MM henceforth) comprise 8 per cent of all remnants 
(■Rhofc Petre..l998.^ . 



Rho & Petre also found that MM SNRs preferen- 
tially occur toward the G alactic midplane where the gas is 
denser . Rho fc Petre cite [ Frail et al.1 (Il99(tl : iReach fc Rhol 
(ll99d) : lKoo fc MoonI (Il997f) to demonstrate that the shocks 
from MM SNRs could be interacting with molecular cloud 
material, suggesting that MM SNRs occur in denser en- 
vironments. In this paper, we examine that hypothesis. 
iKawasaki et al.l (|2005l) have proposed all SNRs pass through 
a composite phase. Their argument is based on over- 
ionisation in the majority of MM SNRs. We examine the 
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hypothesis that all SNRs pass through an MM phase with 
our simulations. 

Two existing models for the formation of MM SNRs 
have figured prominently since the recogniti on of this class 
of rem nant. The first model, proposed by IWhite &; Lond 
(|l99J) . has cool, dense clouds that become engulfed by the 
supernova shock and hot gas bubble. Thermal conduction 
and turbulent mixing transfer heat between the bubble and 
clouds, cooling the interior of the SNR and increasing its 
density as the cold gas in the clouds evaporates and mixes 
with the remnant. 

The second model by ICox et al.l ([l999) proposed that 
thermal conduction acting by itself can cause a SNR to be- 
come centre-bright in x-rays if the remnant is expanding 
into a dense medium. As the sound speed of the thermal 
electrons is significantly greater than the sound speed of the 
ions, the electrons can propagate faster than the outer shock 
and carry energy away from t he interior. The interior density 
increases as the bubble co ols (jSlavin fc Coxlll992 : ICox et al.l 
I1999I : IShelton et al.l ll999V The transport of heat from the 
hot interior of the remnant to the surrounding medium leads 
to lower temperatures and higher densities inside the SNR, 
leading to greater emission in x-rays from the centre of the 
remnant. 

In each of these models, thermal conduction is invoked 
to transport energy from the hot gas to cool gas. They differ 
primarily in how they distribute the cold, dense gas compo- 
nent. In either of these scenarios, the higher densities and 
cooler temperatures of 10^ — 10^ K in the hot gas bubble 
favor the production of soft x-rays in the centres of these 
remnants. We concentrate on the latter model in this work. 

In the presence of a magnetic field, thermal conduction 
becomes anisotropic as the free electrons that are respon- 
sible for the bulk of the heat transport will not propagate 
further than the Larmor rad i us of the ma gnetic field (|SDitzeii 
ll95dlCowie fc McKed[l977l : lBalbusl[l98aV Thus, the geom- 
etry of the magnetic field cannot be neglected when calculat- 
ing the thermal heat flux in a multidimensional calculation. 
As a result, Tillev fc Balsa ra (2006) examined the role of 
anisotropic thermal conduction in the presence of magnetic 
fields and showed that it is important in determining the 
filli ng factor of hot gas in our Galaxy. 

ISlavin fc Coxl (|l992", 1993) used a simplified represen- 
tation of the magnetic fiel d in their early c a lculat ions and 
we improve on that here. I Velazquez et al.l (|2004l) include 
the anisotropy in the thermal conduction introduced by the 
magnetic field, but do not include the Lorentz force of the 
magnetic field on the gas. As a consequence, the magnetic 
field was stirred by the turbulent motions developing behind 
the supernova shock, greatly inhibiting thermal conduction. 
Without the feedback of the magnetic field on the gas, how- 
ever, this can only be regarded as a lower limit on the de- 
gree of thermal c onduction. We removed this restriction in 
iTillev fc Balsara! (|2006), and here we examine further con- 
sequences of removing that restriction. 

A goal of our work, with its limited but interesting set 
of input micro-physics, is to produce diagnostics that may 
be compared with observations. Recent work by the FUSE, 
Chandra and XMM observatories have provided insight 
into the total amount of hot gas in the Galaxy, as traced 
by highly-ionise d stages of oxyge n, esp. O VI and O VI 1 1 
JOegerle 2005: Yao fc Wang 2005: iSavage fc Le hneil l200d) . 



Table 1. Initial conditions of the ISM in our simulations, p is the 
density in amu cni~^; T is the temperature in Kelvin; B is the 
magnetic field in fiG; and L is the grid size in parsecs. 



Run 



B 



LO 


0.7 


8000 


0.0 


75 


LI 


0.7 


8000 


3.0 


75 


MO 


1.0 


10000 


0.0 


50 


Ml 


1.0 


10000 


3.0 


50 


HO 


5.0 


10000 


0.0 


50 


HI 


5.0 


10000 


6.0 


50 



These ions probe gas in the temperature range '-^ 10^ K 
to ^ 10^ K, providing valuable diagnostics of the thermal 
state of the hot material in the Milky Way. While such 
data have so far been used to test global models of the 
ISM (iMcKee fc OstrikeJll977l: iKorpi et al.[ll999l: iKim et al l 
12001'; 'de Avillez & Breitschwerdt' '2004'; 'Balsara et all 
12004; Balsara & Kim 2005; Avillez & Breitschwerdt 2005^ 
Mac Low et all l2005l: iTillev fc Balsaral |2006), it is also 
possible to use high-stage ions as diagnostics for the 
structure of SNRs. Indeed, while earlier observatories 
allowed spectroscopy of small regions of supernova rem- 
nants (e.g., Vedder et al. 1986), recent Chandra and XMM 
observations have provided spatially-resolved m easurements 
of O VII and O viii emission from SNRs (iGaetz et al.l 
| 2000l: iRasmussen et all I2OOII: Ivan der Hevden et al.l l2002l : 
'Flanagan et al.l l2004ri . We further develop the idea of using 
high-stage ions as diagnostics of hot gas in SNRs in this 
paper. 

In this paper we report on the results of a series of 
MHD simulations on the evolution of SNRs with anisotropic 
thermal conduction. We explore a parameter space of initial 
ISM conditions in density, temperature, and magnetic field 
strength (Section [21. In Section |3| we analyse the effects of 
the environment on the evolution of the remnant, in partic- 
ular focusing upon the effects of the environment on pro- 
ducing mixed- morphology remnants. Section 0] explores the 
role that high-stage ions could play in diagnosing the hot 
gas. Section |51 provides some conclusions. 



2 NUMERICAL SETUP 

We carried out our simulations by usin g t he RIEMANN code 
(iRoe fc Balsaral Il996l: iBalsaral Il998al lbi iBalsara fc Spiced 
Il999al lbl: lBalsarall200lL JloS^lEMANN updates the MHD 
equations using a second-order accurate dimensionally un- 
split TVD scheme. We use an implicit solver for the ther- 
mal conduction. We discuss the details of our treatment of 
anisotropic thermal conduction in an acco mpanying paper 
(Balsara et al., in prep a ration). We draw on ISpitzeJ (IIS 
ICowie fc Mc"Ke3 (11977V Balbus ('l986VSla vin fc Coxl E^ 
for our description of the classical thermal conduction flux, 
which is proportional to a temperature-dependent conduc- 
tivity and the projection of the thermal gradient on to the 
local magnetic field: 

Fcias = -aT'^% (b • 




VT 



(1) 

where T is the temperature, a = 6 x 10~^ g cm~^ s~^ K~^^^ 
and b = B/|B| is a unit vector along the magnetic field. 
When the flux becomes saturated, it no longer depends on 
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Figure 1. Simulated x-ray surface brightness in soft and hard x-rays, as described in the text, (a) and (b) show soft and hard x-rays 
from Run LI (low-density ISM) at 20 kyr. (c) and (d) do the same for Run LI at 60 kyr. (e) and (f) show soft and hard x-rays from Run 
HI (high-density ISM) at 20 kyr, while (g) and (h) show the same for Run HI at 60 kyr. The tick marks indicate projected distances of 
10 pc. The colour scale is the logarithm of the surface brightness in ergs s"-*^ cm~^ sr"-*^. 



the value of the temperature gradient, only on the presence 
of a parallel componen t of t he temperature gradient to the 
magnetic field ( Balbu^ lToH^ : 



-50pcfsgn (b • VT) b 



(2) 



whe re Cs is the sound sp eed, p the density and (j) ~ 
0.3 dBalbus fc McKedllQSJ) . A flux limiter is utilized to 
transition from t he classical flux to the saturated flux 
JBalbusfc McKeel [T982: SlaviTXCoi IHi^). The transi- 
tion from the classical (Eq. QJ to saturated (Eq. |2J ther- 
mal conduction causes the flux to ch ange character from 
parabolic to hyperbolic (|Balbuslll986i) . The changing char- 
acter in the flux of thermal conduction is accounted for by 
using a Newton-Krylov multigrid method that correctly in- 
corporates the anisotropic structure of the thermal conduc- 
tion operator in the presence of a magnetic field. We have 
tested our formulation of the thermal con duction against the 
supernova shock radius-time relation of ICioffi et al.l (|l988h 
and found excellent agreement. We also found very good 
agreement with the instability growth rates predicted by 
IPieldl (|l965r l . These tests are described in detail in (Balsara 
et al., in preparation). 

Radiative cooling is calculated using the cooling func- 
tion of MacDonald & Bailev ( 1981). We balance the cooling 
rate of our initial conditions with a constant heating rate, 
such that the initial conditions would remain at constant 
temperature in the absence of the supernova remnant. Our 
choice of equilibrium cooling permits us to do a large num- 
ber of runs without carrying the large number of species that 
would be necessary in a formulation that retains the non- 
equilibrium ionisation (jSlavin fc Coxlll992r i. However, more 
economical formulations for the i nclusion of non-equilib rium 
cooling are becoming available (Beniamin et al.ll2001r ) and 
it is our intent to incorporate such physics in future mod- 



els. The focus of our present paper is to examine the role of 
thermal conduction, which does not depend sensitively on 
the details of the radiative cooling. 

We use a uniform mesh in cylindrical coordinates of 
384x384 zones in the R-z plane to update the MHD equa- 
tions. We chose our initial densities and temperatures of the 
quiescent environment so that they are representative of in- 
teresting interstellar medium conditions. Our lowest density 
runs, labelled "LO" and "LI", have a density of 0.7 amu 
cm~^, a temperature of 8000 K and a magnetic field of 
/iG and 3 /xG, respectively. Our intermediate runs, labelled 
"MO" and "Ml", have a density of 1.0 amu cm~^, a temper- 
ature of 10000 K and a magnetic of /iG and 3 /iG, respec- 
tively. Our highest density runs, labelled "HO" and "HI", 
have a density of 5.0 amu cm~^, a temperature of 10000 K 
and a magnetic field of /xG and 6 /iG, respectively. These 
initial conditions are summarized in Table Q In all of these 
cases, we initialize our supernova with 10^^ ergs of energy, 
partitioned such that one-third of the supernova energy is 
contained in thermal energy, and the other two-thirds in ki- 
netic energy. We initialize our remnant with a radius of 10 
zones in all of our simulations. In order to follow the com- 
plete evolution of the remnant, we wish to ensure that the 
outer shock remains on the grid at all times. The physical 
size of our grid in each of the simulations in Table Q reflects 
this. To serve as a counterpoint we have carried out a com- 
plementary set of simulations without thermal conduction, 
but including radiative cooling. 



3 MIXED MORPHOLOGY REMNANTS 

The structure of supernova remnants that include 
anisotropic thermal conduction in magnetized ISMs has 
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Figure 2. Simulated synchrotron surface brightness oc PB^-^ 
from SNRs, shown at two epochs, (Left) 20 kyr after the explo- 
sion. (Right) 60 kyr after the explosion. The top row illustrates 
the expansion of a SNR into a standard ISM (Run LI). The bot- 
tom row illustrates the expansion of a SNR into a higher-density 
ISM (Run HI). The tick marks indicate projected distances of 10 
pc. Units on the linear colour scale are arbitrary, and scaled such 
that the greatest surface brightness has a value of 256. 
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Figure 3. Radius of the outer shock as it evolves with time, 
showing the difference in the evolution between runs with and 
without conduction, and runs in low or high density ISMs. The 
solid and dashed lines virtually overlie each other. Likewise, the 
dotted and dot-dashed lines also practically overlie each other. 



been studied in lTillev h. Balsaral ^^^. In their fig. 1, they 
intercompare remnants that include thermal conduction to 
those that exclude thermal conduction, but in both cases 
they incorporate the effects of radiative cooling. In the snow- 
plow phase, they find that the temperature of the hottest gas 
is reduced by a factor of '-^ 10 , while the density is increased 
by the same factor, in those remnants that include thermal 
conduction. They also find that the inclusion of a uniform 



magnetic field restricts the thermal conduction in the direc- 
tion of the magnetic field. As a consequence, the structure 
of the hot gas bubble becomes elongated in the direction of 
the magnetic field. 

MM SNRs are primarily diagnosed in x-ray and radio 
emission. Moreover, Chandra has made it possible to image 
the same remnant in soft and hard x-rays with high spatial 
resolution. To examine the influence of thermal conduction 
on the production of the morphology of the remnant at these 
energies, we produce simulated emission maps of our models 
at several epochs during the course of the SNR's evolution. 
We calculate the x- r ay em issivity for 10^ — 10^ K gas using 
iRavmond fc SmithI (|l977l) . These emissivities include the 
contribution from bremsstrahlung, high-order atomic lines, 
recombination and two-photon emission for photons in the 
range of 100 eV - 10 keV. We then integrated the optically 
thin emission along the line-of-sight to create projected in- 
tensity maps of the simulated SNRs in x-rays. 

The projected soft and hard x-ray maps are shown in 
Fig. Q To retain some concordance with Chandra data, we 
define 'soft x-rays' as those with photon energies in the range 
of 300-800 eV, and 'hard x-rays' as those with energies in the 
range of 1-5 keV. The resultant images were not convolved 
with an instrumental spatial or spectral response function. 
Fig. C] (a) and (b) show the soft and hard x-rays in the low- 
density run LI at 20 kyr. Fig. Q (c) and (d) do the same 
at 60 kyr. Fig. Q] (e) through (h) show similar data for the 
high-density run HI at the same epochs. Fig. C](a) and (b) 
show that at early epochs, the low-density run LI is shell- 
bright in soft x-rays, and also shell-bright in hard x-rays. The 
magnetic field draws heat from the hottest central gas and 
carries it to the pole, making Fig. Q^b) brighter at the pole, 
while Fig. QJa) is brighter at the equator. Fig. Q] (e) and (f ) 
show that at 20 kyr the high- density run HI is shell-bright 
in the soft x-rays, bu t is already centre-bright in the hard 
x-rays. IShelton et al.l (|1999.'l also present simulated hard x- 
ray images of remnants expanding into an unmagnetized 
medium of similar density, finding that they too produce 
centre-bright images in the hard x-ray and thereby produc- 
ing a concordance between our data and theirs. 

At later epochs, i.e. 60 kyr. Fig. Q (c) and (d) show 
that the low-density run LI is still shell-bright in soft x-rays, 
but is becoming centre-bright in hard x-rays. Fig.C](g) and 
(h) show that the high-density run HI has become centre- 
bright in both soft and hard x-rays. This shows that thermal 
conduction plays an important role in making model HI 
centre-bright in x-rays. A similar conclusion emerges when 
we cross-compare models LO and HO without magnetic fields. 
The difference between models that include magnetic fields 
and those that exclude it is that without magnetic fields 
thermal conduction operates in both the R and z directions, 
thereby increasing the rate at which energy is transported 
from the hot gas bubble. As a result, thermally conducting 
remnants that exclude magnetic fields make their transition 
to a centre-bright x-ray morphology at an earlier epoch. 

The morphological differences between Fig. (a) and 
(b) require further explanation. The anisotropy in the ther- 
mal conduction flux caused by the presence of the magnetic 
field leads to the incomplete appearance of the shell. Tem- 
perature variations perpendicular to the field are not re- 
duced due to thermal conduction. As the magnetic field is 
swept up by the SNR, the magnetic flux becomes concen- 
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trated in a thin shell. Most of the magnetic field in the inte- 
rior of the remnant nevertheless retains some memory of its 
original orientation. As a result, it is preferentially attached 
to the exterior medium closer to the z-axis. Consequently, 
the outer regions of the hot bubble near the axis of symme- 
try have a lower density and higher temperature than the 
outer parts of the bubble near the midplane, resulting in 
lower emission in soft x-rays and a higher emission in hard 
x-rays near the pole at early epochs. A remnant going off 
into a denser medium is smaller and reaches its snowplow 
phase faster. As a result, thermal conduction is more effi- 
cient at transferring energy within such a remnant. For that 
reason, we see that HI has already become centre-bright 
in hard x-rays, with the soft x-rays showing an incomplete 
shell, see Fig.^e) and (f). We see, therefore, that the choice 
of wavebands used to image the remnant can also make a 
difference in its classification as a centre-bright remnant in 
x-rays. In particular, remnants expanding into a denser ISM 
would be much more likely to be classified as centre-bright 
if they are imaged in hard x-rays instead of soft x-rays. 

Magnetic fields are compressed by the outer shock in an 
SNR, and can even be a mplified if the shock interacts with 
interstellar turbulence (jBalsara et al.l 120011 ^ . Cosmic rays 
can also be swept up by the outer shock and are most likely 
accelerated by it. We do not include cosmic ray acceleration 
processes here. As a result, we calculate synchrotron emissiv- 
ities as Csynch oc PB^^ (seel 
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[19991). The simulated synchrotron maps are illustrated in 
Fig. 12 In both the high-density and low-density environ- 
ments the concentrated magnetic field appears as a bright 
shell behind the shock at all times. Thus, we find a clear 
difference between the morphological appearance at 60 kyr 
in x-rays and radio emission. 

The previous paragraph has shown us that the radio 
morphology is always shell-bright. We have also seen that 
the transition to a centre-bright x-ray morphology depends 
to some extent on the wavebands used. If we use soft x-rays, 
we find that remnants expanding into a denser environment 
enter a mixed-morphology phase within about 60 kyr, while 
remnants expanding into a standard ISM do not transition 
to a mixed morphology within that time. However, even the 
remnants expanding into a standard ISM do transition to 
a mixed morphology after about 90 kyr. This might, how- 
ever, be too late in the life of a remnant for it to be visi- 
ble. We note that the soft x-ray flux has already decreased 
by an order of magnitude between Fig. ^a) and Fig. ^c). 
We expect it to decrease by a similar magnitude by 90 kyr. 
If we use hard x-rays, we find that the remnant expand- 
ing into a denser environment is already centre-bright by 
about 20 kyr, while the remnant expanding into a standard 
ISM does not become centre-bright until about 60 kyr. We 
see, therefore, that the classification as a mixed-morphology 
remnant might depend strongly on the waveband used. We 
can, nevertheless, say with confidence that remnants ex- 
panding into denser ISMs might spend a significant fraction 
of their observable lives as mixed-morphology remnants. We 
cannot say the same for remnants expanding into a stan- 
dard ISM. As a result, we conclude that the conjecture in 
iKawasaki et al.l (|2Q05|) . that all remnants pass through a 
mixed- morphology phase in their evolution, may in theory 
be quite plausible. However, remnants expanding into ISMs 
with standard densities or below-average densities may not 



display themselves as mixed-morphology remnants during a 
significant phase of their observable lifetimes. 

It is also important to realize that thermal conduction 
does not reduce the total energy, it only transports it from 
one place to another. Consequently, our numerical formu- 
lation is also fully conservative. Fig. shows the location 
of the outer shock along the R-axis as a function of time. 
There is clearly no difference in the size of the shock whether 
thermal conduction is included or not. Since the location 
of the outer shock (at least in the radial direction) is un- 
affected by thermal conduction, we can conclude that the 
post-shock pressure should be the same in runs with and 
without conduction. The interior of the hot gas bubble is 
always hot enough (and the sound speed in it high enough) 
to ensure that the bubble achieves pressure balance with 
the post-shock gas. In the presence of thermal conduction 
the hot gas bubble can transport part of its thermal energy 
to some of the swept-up post-shock gas. This process can, 
therefore, lower the temperature in the hot gas bubble. To 
maintain constancy of pressure, the density in the hot gas 
will be correspondingly higher. 



4 HIGH-STAGE IONS AS DIAGNOSTICS OF 
THE HOT GAS 

In this Section, we examine the role of high-stage ions as di- 
agnostics of hot gas in SNRs. Highly-ionised oxygen is read- 
ily observed in X-rays (e.g., through the O vii "triplet" at 
21.60, 21.80, and 22.09 A and O viii Lyman-a at 18.98 A) 
and in the far- ultraviolet (using the O VI doublet at 1033.8 
A). We will calculate the emission from these ions which 
probe a range of temperatures (O VI peaks at '^ 3 x 10^ K, 
O VII at 2.8 X 10^ - 1.8 x 10^ K, and O viii at - 2.2 x 10^ 
K). Thus, these ions provide good diagnostics of gas in these 
temperature regimes, which are important during the rem- 
nant evolution. 

We use the line emissivity calculations fro m the Astro- 
physical Plasma Emission Code (APEC; Smit h et al.ll200lh 
with atomic data from the Astrophysical Plasma Emis- 
sion Database (APED) included in the ATOMDB distribu- 
tion (v. 1.3.1).^ The se calculations assum e collisional ion- 
isation equilibrium (JMazzotta et al.l Il998i) , and we adopt 
an abundance of log 0/H= —3.34 from the solar system 
). In what follows, O VI emission repre- 



(lAsphmd et al.l l2n 
sents the sum of the O VI doublet emission, while O VI i refers 
to the summed emission from the aforementioned triplet; 
O VIII emission refers to Lyman- a 18.98 A emission. We have 
tested the temperature- dependent line emissivities from the 
APEC models using a treatment of O VI emission following 
Is hull fc Slavin fl994j, and we find very good agreement. 

Fig.gjandlSlshow maps of O VI, O vii, and O viii at var- 
ious epochs in simulations of SNRs going through low and 
high density ISMs. The results shown in Fig.^Jinclude ther- 
mal conduction. The results in Fig.|5]do not include thermal 
conduction. In each instance we show simulated images of 
O VI, O VII, and O viii from left to right. Fig. gja) and 
Fig. |5Ja) show these data in the low-density run at 20 kyr, 
and Fig. Elb) and Fig. |5Jb) show the same simulation at 60 



^ Available via _http: / /cxc. harvard. edu/atomdb/| 
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Figure 4. Surface brightness due to emission from oxygen ions 
for (a) Run LI at 20 kyr, (b) Run LI at 60 kyr, (c) Run HI at 20 
kyr, and (d) Run HI at 60 kyr. The ions are (left) O VI, (centre) O 
VII, and (right) O Vlll. The tick marks indicate projected distances 
of 10 pc. Because the distribution of O VI emitting gas originates 
from a very narrow shell, the image shows striations that arise 
from the pixellation of this boundary in the simulations. The 
colour scale is the logarithm of the surface brightness in ergs s~^ 




Figure 5. Surface brighness due to emission from oxygen ions for 
simulations without thermal conduction, (a) Run LI at 20 kyr, 
(b) Run LI at 60 kyr, (c) Run HI at 20 kyr, and (d) Run HI at 
60 kyr. The ions are (left) O VI, (centre) O Vll, and (right) O Vlll. 
The tick marks indicate projected distances of 10 pc. Because the 
distribution of O VI emitting gas originates from a very narrow 
shell, the image shows striations that arise from the pixellation of 
this boundary in the simulations. The colour scale is the logarithm 
of the surface brightness in ergs s"-*^ cm~^ sr"-*^. 



kyr. Similarly, Fig.^Jc) and Fig.|5Jc) show these data in the 
high-density run at 20 kyr, and Fig.^Jd) and Fig.|5Jd) show 
the same simulation at 60 kyr. 

We see that O VI is always concentrated at the bound- 
ary of the hot gas bubble, thereby showing its utility as a 
tracer of the gas that it transitional between the hot and cold 
ph ases, and as a tracer of tur bulent mixing layers as s hown 
bvBorkows ki et al.l(ll99(t l: lBegelman fc McKel(ll99(i) . The 
O VII shows itself to be shell-bright in all instances, with 
and without thermal conduction, indicating that it is not 
a good tracer of the role of thermal conduction in SNRs. 
Physically, this may be due to the broad range of temper- 
atures over which O Vll is the dominant ionisation stage of 
oxygen. By contrast, O VIII probes a narrow range of higher 
temperatures, and the distribution of O Viii is most interest- 
ing. Fig. Eld) shows that the O Vlll is centrally brightened 
at later stages in the evolution of the SNR when thermal 
conduction is included. This trend occurs to a lesser extent 
in Fig. 2Jb), with more O Viii emission arising in the cen- 
tral region when thermal conduction is included. Fig. |5Jb) 
and (d) show a strongly edge-brightened appearance in O 
VIII, even at late epochs, when thermal conduction is not 



included. Thus, thermal conduction plays an important role 
in cooUng down the interior of the hot gas bubble at late 
epochs. This changes the ionisation balance, driving the in- 
terior from O IX to O Viii and making the centre of the 
remnant brighter in O Viii. Comparing Fig. 0] with Fig. |5| 
we see that thermal conduction played a major role in mak- 
ing the SNRs centre bright in O Viii at later epochs in dense 
ISMs. Our simulations therefore demonstrate the utility of 
using O VIII as a diagnostic of remnant age. We hasten to 
point out that our present conclusions for the high-stage 
ions are limited by our use of equilibrium cooling. The ioni- 
sation fractions of the different oxygen ions could change in 
the presence of non-equilibrium physics, and a future round 
of simulations will examine whether this trend is preserved 
when such physics are included. 



5 CONCLUSION 

Anisotropic thermal conduction decreases the average tem- 
perature of hot gas by half an order of magnitude over a 
period of a few million years, and increase the mean den- 
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sity by a similar amount. The temperature in the hottest 
portions of the hot gas bubble is reduced by an entire order 
of magnitude when thermal conduction is included, with a 
corresponding increase in the density. These changes affect 
the emission in high- stage ions rather strongly. In particular, 
the distribution of O viii in a denser ISM transitions from 
shell-bright at early epochs to centre-bright at later epochs 
when the physics of thermal conduction is included. 

In denser ISMs (such as in proximity to giant molecu- 
lar clouds, in the bulges of spiral galaxies and in starburst 
galaxies) the higher density and lower temperatures lead to 
the development of centre-bright soft x-ray emission at late 
stages in the evolution of the SNR. This confirms that it 
is possible to produce MM SNRs in denser environments 
through the process of thermal conduction a lone. Our re- 
sults are concordant with the observations of iRho fc Petrd 
(119981) . who find that the majority of MM SNRs are found 
around dense molecular gas . 

We have also shown that the choice of x-ray waveband 
in which the remnant is imaged can play an important role 
in deducing whether the remnant is in a mixed-morphology 
stage or not. If soft x-rays are used, a remnant expanding 
into a dense ISM will not show itself as centre-bright until 
rather late in its evolution. On the other hand, the same rem- 
nant would show a mixed morphology at reasonably early 
epochs when imaged in hard x-rays. Remnants expanding 
into a less dense ISM do not appear centre-bright in either 
soft or hard x-rays until rather late in their evolution. Even 
at those epochs, the hard x-rays would more likely show a 
centre-bright structure than soft x-rays. 
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